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Abstract 

This paper describes a fully integrated aerody- 
namic/dynamic optimization procedure for helicopter 
rotor blades. The procedure combines performance and 
dynamics analyses with a general purpose optimizer. 
The procedure minimizes a linear combination of power 
required (in hover, forward flight, and maneuver) and 
vibratory hub shear. The design variables include 
pretwist, taper initiation, taper ratio, root chord, blade 
stiffnesses, tuning masses, and tuning mass locations. 
Aerodynamic constraints consist of limits on power 
required in hover, forward flight and maneuver; airfoil 
section stall; drag divergence Mach number; minimum 
tip chord; and trim. Dynamic constraints are on fre- 
quencies, minimum autorotational inertia, and maxi- 
mum blade weight. The procedure is demonstrated for 
two cases. In the first case the objective function 
involves power required (in hover, forward flight, and 
maneuver) and dynamics. The second case involves 
only hover power and dynamics. Hie designs from the 
integrated procedure are compared with designs from a 
sequential optimization approach in which the blade is 
first optimized for performance and then for dynamics. 
In both cases, the integrated approach is superior. 
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Symbols 


AI 

ns 2 

autorotational inertia, X W. r. (lbm-in^) 
j=l 

minimum autorotational inertia (lbm-in^) 

^min 

c d 

airfoil section drag coefficient 

c <k 

maximum allowable section drag coefficient 

V 

Cd max 

largest section drag coefficient at azimuthal 

C D 

angle y 

rotor coefficient of drag 

C L 

rotor coefficient of lift 

c r 

root chord (in) 

c t 

tip chord (in) 

Vn 

minimum tip chord (in) 

DVj 

ith design variable 

DV i 

normalizing factor for ith design variable 

EI xx 

chordwise bending stiffness (lb-ft^) 

EI zz 

flapwise bending stiffness (lb-ft^) 

f i 

ith frequency (per rev) 

f il 

lower bound on ith frequency (per rev) 

f iu 

upper bound on ith frequency (per rev) 

F 

objective function 

®i 

ith constraint function 

GJ 

torsional stiffness (lb-ft^) 

ITER 

number of trim iterations 

n^max 

maximum number of trim iterations allowed 

k i 

ith weighting factor in objective function 

m i 

ith tuning mass (lbm) 



N number of blades 

n integer 

nP frequency, integer multiple of rotor rotational 

frequency (per rev) 

ns number of structural segments 

NDV number of design variables 

P a power available (hp) 

P h power required in hover (hp) 

Phj. f reference power required in hover (hp) 

P ff power required in forward flight (hp) 

P ff f reference power required in forward flight (hp) 

P m power required in maneuver (hp) 

P m reference power required in maneuver (hp) 

Ul ref 

R blade radius from center of rotation (in) 

r distance along blade from center of rotation (in) 

rj distance from center of rotation to center of jth 

segment (in) 

4 per rev nonrotating vertical hub shear in 
forward flight (lbf) 

S 4 reference 4 per rev nonrotating vertical hub 

ref 

shear in forward flight (lbf) 

N per rev vertical nonrotating hub shear in 

forward flight (lbf) 

Svj reference N per rev nonrotating vertical hub 

^ref 

shear in forward flight (lbf) 
tj ith torsional frequency (per rev) 

W total blade weight (lbm) 

Wmax maximum blade weight (lbm) 

Wj total weight of jth structural segment (lbm) 

yj location of ith tuning mass 

yjj. point of taper initiation 

Af increment used in frequency window (per rev) 

9 tw maximum pretwist (deg) 

'P azimuth angle, zero over tail (deg) 

Introduction 

The multidisciplinary nature of the helicopter rotor blade 
design process involves several disciplines including aerody- 
namics, dynamics, structures, and acoustics. Most rotor 
optimization has been applied to single disciplines. For 
example, rotor dynamic optimization is discussed in Refs. 
1-6. Rotor structural optimization is discussed in Refs. 1 
and 7. Rotor blades are designed for optimum performance 
in Refs. 1, 8 and 9. Recently, techniques and strategies for 
integrating disciplines in rotoreraft design procedures have 
been emerging. Such a plan is described in Refs. 10 and 11. 

Progress has been made in developing integrated proce- 
dures for rotor blade design (Refs. 12-14). Reference 12 
describes an optimization procedure which designs a rotor 
blade for combined aerodynamics, dynamics, and structures 


in stages. The blade is first designed for aerodynamic perfo- 
rmance with power required in hover as the objective func- 
tion, chord and twist as design variables, and a constraint on 
the autorotational inertia of the blade. The resulting opti- 
mum design results in a shift in the blade natural frequencies 
due to changes in the blade chord and twist distributions. A 
blade-frequency placement optimization is then performed to 
bring the blade natural frequencies to within prescribed 
windows. Finally, the rotor power is minimized with all 
the constraints (frequency, stress, fatigue life, and aeroelastic 
stability). References 13 and 14 describe the formulation of 
a multidisciplinary approach to rotor blade design for impr- 
oved performance and reduced fuselage vibrations. The objec- 
tive function is a linear combination of power required in 
hover, power required in forward flight, and vibratory load. 
The design variables are perturbations from the initial design 
of the following: linear and/or nonlinear twist and chord 
distributions, stiffnesses (flapwise, chordwise, and tor- 
sional), section mass and section moment of inertia. The 
constraints are a lower bound on blade weight for minimum 
autorotational inertia capability, an upper bound on section 
angle of attack at the blade tip, and a lower bound on tip 
chord. The rotor solidity and frequencies are constrained to 
be close to the initial values. The constraint on solidity is 
imposed so that the optimized design will have tire same 
performance as the initial blade without compromising 
maneuvering capability. Reference 14 extended this work to 
included additional design variables - blade sweep and offsets 
of the center of gravity from the aerodynamic center. 

In the present work, a fully integrated aerody- 
namic/dynamic optimization procedure is described. The 
procedure accounts for aerodynamic performance and 
dynamics simultaneously. In this work, constraints are 
enforced on power required, airfoil section stall, drag diver- 
gence Mach number, minimum tip chord, trim, rotor frequ- 
encies, autorotational inertia, and blade weight. This pro- 
cedure accounts for the interactions between disciplines by 
simultaneously changing the design variables to satisfy the 
design requirements and optimizing a composite measure of 
both performance and dynamics. The procedure is similar to 
that of Refs. 13 and 14 but differs in several ways. First, in 
the present work, stiffness distributions, tuning masses, and 
tuning mass locations are design variables. Second, a man- 
euver flight condition is included in the objective function 
and constraints. Third, no constraint on solidity is neces- 
sary since the the blade is designed for a constant lift in 
forward flight and a constant lift in maneuver. Fourth, the 
frequency constraints are formulated so that the frequencies 
are away from integer multiples of the rotor speed and are 
not confined to be near the initial frequency values. 

This paper describes the application of the optimization 
procedure to the design of a l/6th scale model of a rotor 
blade for a utility helicopter. Results for two composite 
objective functions are presented - one involving perfor- 
mance (in hover, forward flight, and maneuver) and 
dynamics; the oilier involving performance in hover and 
dynamics. For both cases comparisons are made between 
integrated and sequential optimization approaches to 
rotor blade design. The integrated approach accounts for the 
interactions between disciplines by simultaneously changing 
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both aerodynamic and dynamic design variables and includ- 
ing both performance and dynamics in the objective function 
and constraints. The sequential approach is to optimize the 
blade first for performance using only aerodynamic design 
variables and performance constraints. This performance- 
optimized design is then optimized for dynamics by adjust- 
ing the blade stiffnesses and adding tuning masses. The 
paper demonstrates that the integrated optimization approach 
is generally better than the sequential approach from the 
standpoint of better designs. 


Rotor Blade Design Considerations 

In the present work the rotor blade is designed for per- 
formance and dynamics for three flight conditions: hover, 
forward flight, and maneuver. The maneuver flight con- 
dition simulates a sustained pull-up maneuver in terms of a 
load factor on the forward flight lift requirement. 
Satisfactory aerodynamic performance is defined by the fol- 
lowing four requirements. First, the power required for any 
flight condition must be less than the available power. 
Second, airfoil section stall along the retreating side of the 
rotor disk must be avoided and the section drag divergence 
Mach number on the advancing side of the rotor disc must 
not be exceeded. The stall requirements are handled by requi- 
ring the airfoil sections distributed along the rotor blade 
operate at section drag coefficients less than a specified 
value, . Tltircf the rotor must trim at each flight condi- 
tion. Fourth, the blade tip chord must be larger than a 
prescribed minimum value, c. 

min 

For this work, satisfactory dynamics is defined in terms 
of limits on vibrational frequencies. The blade is designed 
so that the natural frequencies (both bending and torsional) 
are away from integer multiples of the rotor speed. The 
blade must have sufficient autorotational inertia in case of 
engine failure. In addition to satisfying these design requir- 
ements, the blade weight must not exceed some upper limit, 

^max- 


Optimization Formulation 
Design Variables 

The nineteen design variables shown in Fig. 1 consist of 
aerodynamic quantities describing the blade planform and 
dynamic quantities describing the blade structural properties. 
The four aerodynamic design variables are the point of taper 
initiation y^, root chord c r , taper ratio Cj/c t , and maximum 

pretwist 0 tw . The blade is rectangular to y^ and then tapers 

linearly to the tip. The pretwist (blade structural and aero- 
dynamic twist are assumed to be the same) varies linearly 
from the center of rotation to the tip. Nine dynamic design 
variables include the blade chordwise, flapwise, and torsional 
stiffnesses (denoted by EI XX , EI ZZ , and GJ, respectively) at 
the blade root, point of taper initiation, and blade tip. The 
stiffnesses are assumed to vary linearly between these three 


points. Although the three stiffnesses are treated as inde- 
pendent variables, they are in fact related at a given cross 
section by a single set of cross-sectional dimensions and 
material properties. However, in the present formulation 
this relationship is not incorporated. This reconciliation of 
stiffnesses and their physical realizability are assumed to be 
relegated to a separate design level where the cross section of 
the blade is determined. In principle, this multilevel 
decomposition of the rotor blade design is described in Ref. 
11. A demonstration of the feasibility of this type of 
decomposition where independent global design variables are 
determined at one level and the reconciliation of the relation- 
ships among the global variables using detailed design 
variables is done at a lower level is described in Ref. 15. In 
the absence of the lower level the mutual dependence of the 
stiffnesses and physical realizability of the corresponding 
cross sections is accounted for by initialization to an exist- 
ing blade and the expectation that these variables will not 
depart much from their initial values. The remaining six 
dynamic design variables are three tuning masses (denoted 
by mj, m 2 , m 3 ) and their locations (denoted by yj, y 2 , and 
y 3 ). The total blade mass consists of the structural mass 
(which is assumed constant) plus the sum of the tuning 
masses. It is assumed that the center of gravity and aerody- 
namic offsets are coincident with the blade elastic axis. The 
number of blades, rotor radius, rotational velocity, airfoils, 
and airfoil distribution are preselected and fixed. Upper and 
lower bounds on the design variables are summarized in 
Table 1. 

Constraints 

In this section of die paper, the rotor blade design requir- 
ements are expressed as mathematical constraints. By con- 
vention a constraint function denoted by gj is satisfied if it 
is less than or equal to zero. The constraints are grouped 
into performance constraints and dynamic constraints. The 
performance constraints are imposed for all three flight con- 
ditions. The dynamic constraints are imposed only in for- 
ward flight and maneuver. Parameters used in the con- 
straints are summarized in Table 1 . 

Performance Constraints - Recall die performance 
constraints are on power required, trim, stall, and blade tip 
chord. The requirement that die powers required in hover, 
forward flight, and maneuver be less than the power avail- 
able translates into three constraints 

gj= P f /P a - 1 < 0 for each flight condition ( 1 ) 

where P r and P a are the power required and die power avail- 
able, respectively. 

The requirement dial die airfoil sections not stall in for- 
ward flight and maneuver and dial the drag divergence Mach 
number be avoided translates into constraints on die airfoil 
section drag coefficient, c^. This leads to 24 constraints per 
flight condition since the c^'s are evaluated at every 15 

degrees around die azimuth. At a given azimuthal angle 'P 
die constraint is formulated as 
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( 6 ) 


g r = C J /c d - 1.S0 'P=15,30,45,...,360 (2) 

u max u all 

'P 

where c,, is the allowable drag coefficient and c, is the 
all "max 

largest drag coefficient along the blade radius outside the 
reverse flow region at the azimuthal angle V P. 

In this work, an isolated rotor analysis is used which 
trims the rotor to constant lift Cj^ and drag Cq and zero 
flapping angle relative to the shaft using collective, lateral 
cyclic and longitudinal cyclic pitch. Trimming to a con- 
stant lift ensures that the rotor has no loss in lift capability 
even if solidity decreases. The trim requirement is difficult 
to translate into a mathematical constraint. The trim con- 
straint in forward flight and maneuver is implemented using 
the method developed in Ref. 8 which expresses the con- 
straint in terms of the number of trim iterations ITER, the 
maximum number of trim iterations allowed riER max> and 

the design variables DVj. The heuristic trim constraint is 
given by 


NDV 


gi = (lTER-ITER max+ l)( X ) S 0 (3) 

P=1 % 


where NDV is the number of design variables, DVp is die 

* 

pth design variable and DVp is a normalizing factor so that 

the each term in the summation is approximately equal to 
one and positive. In development of Eqn. 3 in Ref. 8, it 
was found that the addition of the summation term helped to 
improve convergence by imparting additional information to 
use in the search direction for new design variable values. 


The final performance requirement is a constraint used to 
ensure that the blade tip chord does not become too small 


gi = 1 - c t /c t . 

1 l min 


(4) 


where c. is the tip chord and c, is the minimum tip 

l min 

chord allowed. These four design requirements were used in 
die performance optimization procedure of Refs. 8 and 9 and 
are considered die performance constraints. 

Dynamic Constraints - Recall the dynamic 
constraints are on frequencies, total blade weight, and 
autorotational inertia. The constraint on the kth frequency 
fj, (either a bending or a torsional frequency) is formulated 
such tiiat the frequency is away from integer multiples of 
the rotor speed by an amount Af 

f k 

gj = “■ - 1 <0 (5) 

leu 
and 


f k 

gi = 1 - < 0 
lei 

where f has a value that is Af below n+1 per rev and 

has a value that is Af above n per rev for the applicable n. 
For example, suppose Af is 0.1 per rev and is 5 6 per 
rev, then nP would be 5 per rev and (n+l)P would be 6 per 
rev. Thus f 4u and f^j would be 5.9 per rev and 5.1 per rev, 

respectively. Formulating the constraints in this manner 
allow the frequencies to change provided they avoid 
approaching integer multiples of the rotor speed. This 
formulation is different from the approaches used in Refs. 
12-14 where the frequencies are kept within prescribed win- 
dows close to the reference blade frequencies. In this work 
constraints are placed on frequencies in both forward flight 
and maneuver. Although these are the natural frequencies 
and should be independent of flight condition, the blade col- 
lective pitch may be different in forward flight and maneu- 
ver. Thus the amount of modal coupling would vary and it 
is possible that the frequencies could be different for the two 
flight conditions. 

The constraint that the blade weight be less than some 
maximum value is formulated as follows 

gi = W/Wmax -ISO (7) 

where W is the total blade weight and W max is the maxi- 
mum allowable weight. The total blade mass is determined 
from the structural mass distribution (which is assumed 
constant) plus the sum of the tuning masses. 

Finally, the blade must have enough autorotational iner- 
tia AI for safe autorotation in case of engine failure. The 
constraint is formulated so that the autorotational inertia of 
the blade is greater than some minimum value AI m - n 

gi = 1 - AI/AI min < 0 (8) 


Objective function - The objective function used in 
this work is similar to the objective function defined in 
Refs. 10 and 11 which was a linear combination of power 
required (in hover, climb, forward flight, and maneuver) and 
hub shear. In the present work the climb term has been 
omitted. The objective function to be minimized is a com- 
bination of performance (the power required for each flight 
condition) and dynamics measure (the N per rev nonrotating 
vertical hub shear in forward flight where N is the number 
of blades) and is formulated as follows 


Ph 

F = k, tt- 0 - + fc 


^ P, 


a ref 


2 P, 


p ff 


+ k 


ff. 


ref 


3 P, 


m 




•+ k 


m 


ref 


4 S 


ff 


N, 


(9) 


ref 


where Pj r P^-, and P m are the powers required in hover, 
forward flight, and maneuver, respectively. Sj^ is the N 
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per rev nonrotating vertical hub shear in forward flight. The 
terms ki, k^, k?, and k 4 are weighting factors. P h , 
J ^ "ref 


p ff „> P m j and Sxj are reference values used to normal- 
n ref “Vef iN ref 

ize and nondimensionalize the objective function compo- 
nents. 


Kotor Analyses 

The analyses used in this work are the Langley-developed 
hover analysis program HOVT (a strip theory momentum 
analysis based on Ref. 16) and the comprehensive helicopter 
analysis program CAMRAD/JA (Ref. 17) for forward flight 
and maneuver. HOVT is used to predict power required in 
hover using nonuniform inflow (no wake is included). 
CAMRAD/JA is used to calculate rotor performance, loads, 
and frequencies. In this work the CAMRAD/JA analyses 
are performed with uniform inflow with empirical inflow 
correction factors. Both HOVT and CAMRAD/JA use 
tables of experimental two-dimensional airfoil data. 


Optimization Methods 

The optimization methods used in this work are the gen- 
eral purpose optimization program CONMIN (Ref. 18) and 
an approximate analysis used to reduce the number of 
HOVT and CAMRAD/JA analyses during the iteration pro- 
cess. CONMIN is a general purpose optimization program 
which uses the method of feasible directions for constrained 
function minimization. The approximate analysis is Used to 
extrapolate the objective function and constraints with linear 
Taylor Series expansions using derivatives of the objective 
function and constraints with respect to the design variables. 
The assumption of linearity is valid over a suitably small 
change in the design variable values and will not introduce a 
large error into the analysis provided the changes are small. 
Errors which may be introduced by use of the approximate 
analysis are controlled by imposing “move limits” on each 
design variable during the iteration process. A move limit 
which is specified as a fractional change of each design vari- 
able value is imposed as an upper and lower design variable 
bound. 


Implementation of Optimization Procedure 

The optimization procedure (Fig. 2) consists of an outer 
loop denoted by “Cycle” and an inner loop denoted by 
“Iteration”. First, preassigned parameters such as the blade 
radius, airfoil distribution, and number of blades are set. An 
optimization cycle is initiated. The aerodynamic and struc- 
tural properties such as twist and chord distributions, radial 
station locations, solidity, blade weight, and autorotational 
inertia are calculated using the current design variable values 
in the box labelled “Design variable preprocessors”. The 
HOVT analysis is then performed to obtain the power requi- 
red in hover. Two CAMRAD/JA analyses (forward flight 
and maneuver) are then performed to obtain the power requi- 
red, trim information, c^'s for the stall constraints, natural 


frequencies, and hub shears. This information is then used 
to formulate the objective function and constraints. Since 
CONMIN and the approximate analysis need derivatives of 
the objective function and constraints, a sensitivity analysis 
is performed to obtain finite difference derivatives of the 
objective function and constraints with respect to the design 
variables. These derivatives are obtained by perturbing each 
design variable one at a time and going through the design 
variable preprocessor, HOVT, and CAMRAD/JA analyses. 
The inner loop consists of CONMIN and the approximate 
analysis. New values for the design variables are obtained 
and the outer loop is re-entered. Convergence is obtained if 
the objective functions from three consecutive cycles are the 

same within a tolerance of 0.5 x 10"^. 


Test Problem 


Shown in Fig. 3 is a l/ 6 th-scale wind tunnel model of a 
rotor blade for a four-bladed utility helicopter. This blade 
has a rectangular planform to 0.80R (80 percent radius) and 
then tapers to the tip with a 3-to-l taper ratio. The blade 
has a radius of 56.22 in and a root chord of 5.40 in. Three 
sets of advanced airfoils are used along the blade. The 
RC(4)-10 airfoil (Ref. 19) is used to 82.5 percent radius, the 
RC(3)-10 (Ref. 20) airfoil is used from 87.5 to 92.5 percent 
radius, and the RC(3)-08 (Ref. 20) airfoil is used from 97.5 
percent radius to the tip. Details of the blade can be found 
in Ref. 21. This blade will be referred to as the reference 
blade . 

An analytical model of the reference blade with 19 aero- 
dynamic segments for HOVT, 50 structural segments and 18 
aerodynamic segments for CAMRAD/JA is used to obtain 
values for the maximum allowable blade weight, the mini- 
mum value for the autorotational inertia, and the structural 
mass distribution. The flight conditions are a constant lift 
of 1-g (331 pounds, Cl= 0.0081), propulsive force of 32 
pounds (C d =- 0.000811), and an advance ratio of 0.35 for 

the forward flight condition and a constant lift of 401 
pounds (Cj_=0 .00985), a propulsive force of 23 pounds 

(Cp)=-0.000596), and an advance ratio of 0.3 for the maneu- 
ver flight condition. The maneuver flight condition is for a 
load factor of 1.22. From the modal analyses in 
CAMRAD/JA using ten bending modes and five torsional 
modes, it is found that only the first six bending frequencies 
are below 10 per rev and need to be constrained for a four- 
bladed rotor. Since fj corresponds to a rigid body mode and 
f 2 is the 1 per rev, the first two frequencies are not con- 
strained. Constraints are placed on the first four bending 
frequencies (fg, £4 and fg flapping-dominated and f 5 lead-lag 
dominated) and die first two torsional frequencies (tj repre- 
senting the rigid body torsional mode due to the control sys- 
tem stiffness and t 2 representing die first elastic torsional 
mode). From the procedure, this blade has a total weight of 
3.05 lbs and an autorotational inertia value of 3411 lbm- 
in^. In this work, the blade is to be designed so that the 
weight is not increased by more dian 15 percent (W max =3.5 
lbs) and the autorotational inertia is increased by at least 1 
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percent (AI m j n =3456 lbm-in^) from that of the reference 

blade. The values for minimum tip chord (c f ), power 

hnin 

available (P a ), and maximum allowable drag coefficient 

(c d ) are 1 in, 20 hp, and 0.1, respectively. A Af of 0.1 
an 

per rev is used for the frequency constraints. Since a four- 
bladed rotor is used as the test problem, the 4 per rev nonro- 
tating hub shear is used for Sjq in the objective function 
given by Eqn 9. Parameters and flight conditions are sum- 
marized in Table 1. 

The initial blade design (the starting point for the opti- 
mization) used in this work is shown in Fig. 4. This blade 
has a rectangular planform with a maximum pretwist of -9.0 
deg and blade root chord of 5.40 in. This blade has the same 
root chord, mass distribution, and stiffness distributions at 
the root, 0.8R (point of taper initiation of the reference 
blade), and tip as the reference blade. The stiffnesses are as- 
sumed to vary linearly between these points. The nonstruc- 
tural mass distribution only depends on the tuning masses 
and their locations. Note that the initial blade does not sat- 
isfy the minimum autorotational inertia requirement. 

Results 



A total of 19 design variables and 80 constraints are 
used. There are 4 aerodynamic design variables ( yy, c f , 
Cj/c^, and 0 tw ) and 15 dynamic design variables (EI XX EI ZZ , 
and GJ at the root, point of taper initiation, and tip; mj; 
m 2 ; m 3 ; yj, y 2 , and y$). The 54 performance constraints 
are: 3 constraints on power required (one per flight condi- 
tion), 2 constraints on trim (one for forward flight and one 
for maneuver), 48 constraints on c^ (24 for forward flight 
and 24 for maneuver), and one constraint on minimum tip 
chord. The 26 dynamic constraints are: 24 frequency con- 
straints (12 for forward flight and 12 for maneuver), a con- 
straint on total blade weight, and a constraint on minimum 
autorotational inertia. 


Case 1 Objective Function - The objective func- 
tion is a combination of the power required (in hover, for- 
ward flight, and maneuver) and the 4 per rev nonrotating ver- 
tical hub shear in forward flight 



The optimization procedure described previously is ap- 
plied to two versions of the objective function given in Eqn. 
9. For each case the normalizing factors P h , P ff , and 

ref ref 

P m are each chosen to be 1 5 hp and S w is chosen to be 

Vef 1 N ref 

2 lbf (based on analysis of the initial blade). The first case 
has an objective function which is a linear combination of 
the power required in hover, forward flight, and maneuver 
and the 4 per rev nonrotating vertical hub shear in forward 
flight. The reference blade was originally designed for per- 
formance. Therefore, the objective function is chosen to be 
one dominated by performance with little emphasis on 
dynamics. Of the three powers it is assumed that it is most 
important to reduce the power required in hover - it will 
have twice the weight as the other two powers. Several 
values were tried for die weighting factor on the hub shear 
term. It was found that to obtain the proper balance be- 
tween performance and dynamics, kq has to be between one 
and two orders of magnitude less than kj. Thus, for this 
case, the weighting factors are chosen to be kj —15.0, 
k 2 =k 3 = 7 . 5 , and kq = 0.025. This objective function will be 
referred to as Case 1 . The second case has an objective func- 
tion which is a linear combination of the power required in 
hover and the hub shear in forward flight (kj=15.0, k 4 = 2 . 0 , 
k 2 =k 3 = 0 . 0 ). There are two reasons for investigating the 
second case. First, including only hover power in the ob- 
jective function is similar to the conventional performance 
optimization described in Refs. 8 and 9. Second, it is of 
interest to see the effect of a larger emphasis on hub shear in 
the objective function. This objective function will be re- 
fereed to as Case 2. 


where Pjj ^ Pff ^ and ^ are each 1 5 hp, and S 4 ^ is 2 

lbf. Results for the initial and integrated-optimized blades 
are summarized in Table 2 . The optimized design has more 
pretwist (-16 degrees) and less root chord (4.38 in) than the 
initial design. The planform has changed from rectangular 
to tapered at 0.68R with a taper ratio of 1.79. The tuning 
masses are located about 0.395R to 0.475R with the two 
largest masses concentrated near 0.40R. The autorotational 
inertia increased 4.4 percent over the initial design 
(representing a 3 percent increase over the minimum re- 
quirement). There is a 9.5 percent increase in blade weight 
(part of the increase is due to the autorotational inertia re- 
quirement). All frequencies are away from per rev values. 
The powers required in hover, forward flight, and maneuver 
are reduced by 6.5, 6.7 and 4.1 percent, respectively, from 
the initial design. The forward flight hub shear is reduced 
from 2.19 to 1.89 lbf (representing a 13.7 percent reduc- 
tion). This is a large reduction in spite of the small empha- 
sis on hub shear in the objective function. The maneuver 
hub shear is reduced slightly from 0.95 to 0.92. Recall 
there is no constraint on maneuver hub shear and only for- 
ward flight hub shear is in the objective function. 

The stiffness distributions for the initial and final designs 
are shown in Fig. 5. As shown in Fig. 5a, inboard of 0.5R 
the chordwise bending stiffness EI XX is much higher than 
the initial design. From this point to about 0.9R, the stiff- 
ness is smaller than the initial design with the smallest 
value around 0.68R (the point of taper initiation). When 
point of taper initiation moved the stiffness distribution 
changed with it. Outboard of 0.9R the stiffness is slightly 
higher than the initial design. Inboard of 0.5R the flapwise 
bending stiffness EI ZZ distribution is higher and outboard of 
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this location it is less than the initial design (see Fig. 5b). 
The torsional stiffness GJ is higher than the initial design 
inboard of 0.7 1R. GJ is smaller until 0.95R when it is 
slightly higher than the initial design (see Fig. 5c). 

Comparison with Reference Blade - It is of in- 
terest to see how this integrated design compares with the 
reference blade for these same flight conditions. As shown 
in Table 3, the analysis predictions for the integrated design 
and the reference blade are similar. The integrated design has 
the same pretwist as the reference blade. Both planforms are 
similar with the integrated design having less solidity than 
the reference blade. The difference in planforms is primarily 
due to the choice of flight conditions. The reference blade 
was designed by nonoptimization techniques for slightly dif- 
ferent flight conditions and design requirements. 

Case 2 Objective Function - In Case 2, only 
hover performance and forward flight dynamics are included 
in the objective function 


F = 15.0 


+ 2 . 0 ; 


+ ff 


( 11 ) 


1 ref 


+ ref 


Integrated versus Sequential Optimization 
Results 


It is of interest to compare the integrated optimization 
approach and a sequential optimization approach. The two 
approaches are summarized in Table 5 using the objective 
functions defined in Eqns. 10-15. The integrated approach is 
a one-step optimization procedure and uses all the aerody- 
namic and dynamic design variables, all the aerodynamic and 
dynamic constraints, and has an objective function which is 
a composite measure of performance and dynamics. The 
sequential approach is a two step optimization procedure. 
The first step consists of a performance optimization similar 
to that of Refs. 8 and 9. The blade is designed using only 
the four aerodynamic design variables (©tw* ytr> c^t* ^ d ® 1 
the performance constraints (power required, trim, stall, and 
minimum tip chord). In Step 2, the optimized-aerodynamic 
design variables are then held constant and the design is 
optimized for dynamics using the only the dynamic design 
variables (stiffnesses, tuning masses, and tuning mass loca- 
tions) and dynamic constraints (frequencies, blade weight, 
and autorotational inertia) The trim constraint is still pre- 
sent . A sequential optimization procedure is done for both 
Case 1 and Case 2. Results are compared with those 
obtained using the integrated design approach. 


where P h is 15 hp, and S 4 is 2 lbf. Results are pre- 
ref ref 

sented in Table 4 for the initial and die optimized designs. 
The optimized design has less pretwist (-8.1 degrees) and a 
smaller root chord (4.5 in) than the initial design. The plan- 
form has changed from a rectangular blade to one that tapers 
at 0.8R with a taper ratio of 1 .36. The tuning masses are 
located between 0.25R and 0.41R with the largest mass at 
0.29R and a slightly smaller mass at 0.25R. There is a 10 
percent increase in blade weight. There is improvement in 
the power required for all flight conditions although not as 
much as in Case 1. The Case 2 design reduces the power 
required in hover and in forward flight by 3.9 and 4.4 per- 
cent, respectively. The two torsional frequencies (tj and t 2 ) 
are at their lower and upper bounds, respectively. The for- 
ward flight hub shear is reduced by 67.6 percent relative to 
the initial design (from 2.19 to 0.71 lbf). The maneuver 
hub shear is increased by 65.3 percent (recall there is no 
requirement on maneuver hub shear either in the objective 
function or constraints). The larger contribution of the hub 
shear term in die objective function seems to have a large 
influence on die resulting design. In this case, there is little 
change in the twist and less taper compared to Case 1. 

The stiffness distributions are shown in Fig. 6 . As 
shown in Fig. 6 a, the chordwise stiffness EI XX is slightly 
higher than the initial design inboard of the point of taper 
initiation and dien decreases to a smaller value at die tip. 
The flapwise stiffness EI ZZ is larger tiian die initial stiffness 
over the blade span (see Fig. 6 b). The torsional stiffness GJ 
is also larger over the entire span with the greatest difference 
at the point of taper initiation (see Fig. 6 c). 


Case 1 Objective Function - In the sequential 
optimization approach, the blade is optimized first for per- 
formance using only the aerodynamic design variables and 
constraints. In Step 1 the objective function in Eqn. 10 
becomes 

Pi, Pff P m 

F = 15.0 -^ JL - + 7.5 + 7.5 (12) 

^ref ^ref m ref 

where Pt, , P ff , and P m are each 15 hp. Theresult- 
u ret ref 11 ref 

ing blade is shown in the column labelled “Performance” in 
Table 6 . The performance-optimized blade has a maximum 
pretwist of -16.0 degrees, is rectangular out to 0.51R, it has 
a taper ratio of 1.65, and a root chord of 4.73 in. This 
design requires 14.4 hp, 12.7 hp, and 1 1.8 hp in hover, for- 
ward flight, and maneuver, respectively. 

Now, this performance-optimized blade is examined from 
the dynamics point-of-view by looking at die blade frequen- 
cies, hub shear, autorotational inertia, and blade weight. 
The blade frequencies have shifted from their original values. 
The second bending frequency (4 (5.04 per rev) is unaccept- 
able (recall frequency windows of 0.1 per rev were used, so 
this frequency should be between 5.1 and 5.9 per rev). Also 
die performance-optimized design has 4 per rev hub shears 
which increased by 32.4 percent (from 2. 19 to 2.90 lbf) and 
52.6 percent (from 0.950 to 1.45 lbf) over the initial values 
for forward flight and maneuver, respectively. The blade 
does not meet die minimum autorotational inertia require- 
ment (recall die initial blade does not meet the minimum 
autorotational inertia requirement). 
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In Step 2, the performance-optimized design is optimized 
for dynamics by minimizing the 4 per rev hub shear in for- 
ward flight using dynamic design variables and dynamic 
constraints with performance design variables held constant. 
The objective function from Eqn. 10 is 


F = 



(13) 


Case 2 Objective Function - To see if the advan- 
tage of the integrated approach is dependent on the choice of 
objective function, the integrated and sequential approaches 
are compared for the case where only hover performance and 
dynamics are included in the objective function. In Step 1, 
the objective function in Eqn. 1 1 becomes 

p h 

p _ _J! (i 4) 


where S 4 is 2 lbf. Results for the dynamic optimization 
ref 

are shown in the column labelled “Dynamics” in Table 6. 
This is the final sequential optimum design. The unaccept- 
able frequency f 4 is brought into the acceptable frequency 
window (5.10 per rev) with a slight shift in the other bend- 
ing frequencies. This sequential optimum design has larger 
4 per rev hub shears than the initial design (by 14.2 percent 
- from 2.19 to 2.50 lbf for forward flight - and by 38.9 
percent - from 0.95 to 1 .32 lbf for maneuver). 

Table 6 also contains the integrated results (from Table 
2). Both the sequential and integrated optimized designs 
have the same pretwist. The integrated approach design has 
the point of taper initiation further outboard and a larger 
taper ratio than the sequential approach design. Both 
designs have the largest tuning mass at about the same loca- 
tion (0.40R). The integrated design has larger tuning 
masses and thus weighs more but has a better autorotational 
inertia capability than the sequential design. Both designs 
reduce power required in hover by the same amount. The 
integrated design requires slightly less forward flight power, 
but more maneuver power. Both designs have acceptable 
frequencies. The integrated design has frequencies in the 
same windows as the initial design. The sequential design 
has two frequencies which have shifted windows (f$ and t 2 ). 
The largest difference between the two designs is in the hub 
shear values. The sequential design increases both forward 
flight and maneuver hub shears by 14.2 and 38.9 percent 
over the initial design values, respectively. The integrated 
design reduces forward flight hub shear by 13.7 percent and 
maneuver hub shear by 3.2 percent over the initial design. 


where Pu is 15 hp. The performance-optimized design 
ref 

shown in the column labelled “Performance” of Table 7 
changed the pretwist from -9.0 to -16.0 degrees, the plan- 
form from rectangular to tapered with the point of taper ini- 
tiation at 0.63R, a taper ratio of 1.28, and a root chord of 
4,04 in. This design required 14.4 hp, 12.6 hp, and 1 1 .9 hp 
in hover, forward flight, and maneuver, respectively. As 
before in the performance optimization step in Case 1, there 
is a shift in blade frequencies resulting in two frequencies 
being too close to a per rev value (fg=8.99 per rev which is 
too close to 9 per rev and 12=14.94 per rev which is too 


close to 15 per rev). As before the autorotational inertia is 
inadequate since the initial blade has inadequate autorota- 
tional inertia. 


In the dynamic optimization of Step 2, the objective 
function in Eqn. 1 1 becomes 


F = 


+ ff 


+ ref 


(15) 


where S 4 is 2 lbf. These results are presented in the col- 
ref 

umn labelled “Dynamics” of Table 7. This final sequen- 
tially-optimized design has satisfactory frequencies. Tuning 
masses are placed between 0.29R and 0.4 1R. Recall in 
Case 1 the tuning masses are more localized. The blade 
weight increases by 3.9 percent from die initial design. The 
autorotational inertia is acceptable. 


The blade stiffnesses are presented in Fig. 7 for the ini- 
tial, sequentially-optimized, and integrated-optimized 
designs. For the sequentially-optimized design, the chord- 
wise bending stiffness EI XX is higher inboard of 0.5R than 
the initial design and less stiff outboard (see Fig. 7a). The 
integrated design EI XX is higher than die sequential design to 
0.45R. The sequential design flapwise bending stiffness 
EI ZZ is slighdy higher inboard of 0.30 R and less outboard 
than the initial design, while the integrated design has an 
EI ZZ distribution higher than die initial distribution until 
around 0.55R (see Fig. 7b). Relative to the initial design, 
the torsional stiffness GJ for the sequentially-optimized 
design inboard of die point of taper initiation is die same, 
while outboard die stiffness has decreased (see Fig. 7c). For 
die integrated design the GJ distribution is higher dian diat 
of die sequential design. 


Integrated optimum results (from Table 4) are included in 
Table 7 in the column labelled “Integrated” for comparison. 
In this case, diere is more of a difference between the 
sequential and integrated designs than in Case 1. The 
sequential design has more twist, a point of taper initiation 
more inboard, less taper, and a smaller root chord dian the 
integrated design. The emphasis on the 4 per rev forward 
flight hub shear in the integrated objective function results 
in less twist and taper. Using die sequential approach, the 
forward flight hub shear is reduced by 20.1 percent over the 
initial design (from 2.19 to 1.75 lbf). The integrated 
approach gives a design with slightly worse performance 
than the sequential design but has reduced die 4 per rev hub 
shear in forward flight 67.6 percent over the initial blade 
compared with a 20.1 percent reduction for die sequential 
approach design. Bodi designs increased die maneuver hub 
shear significantly. 
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As shown in Fig. 8a, the sequential approach chordwise 
bending stiffness EI XX is larger than the initial design stiff- 
ness inboard of 0.45R and outboard of 0.82R. The inte- 
grated approach design chordwise bending stiffness EI XX is 
slightly larger than the initial design stiffness inboard of 

0.9R. In the sequential approach, the flapwise stiffness EI ZZ 
is smaller than the initial blade stiffness, while in the inte- 
grated approach the opposite is true (see Fig. 8b). In the 
sequential approach, the torsional stiffness GJ is the same 
up to the point of taper initiation and then decreases to the 
tip compared to the initial GJ, while in the integrated 
approach the GJ is higher except at the tip (see Fig. 8c). 


Concluding Remarks 

A fully integrated aerodynamic/dynamic optimization 
procedure has been developed for helicopter rotor blades. The 
procedure combines performance and dynamics analyses for 
hover, forward flight, and maneuver with a general purpose 
optimizer. The procedure minimizes an objective function 
which is a composite measure of performance and dynamics. 
Specifically, the objective function is a linear combination 
of power required (for hover, forward flight, and maneuver) 
and vibratory hub shear. The design variables include 
pretwist, taper initiation, taper ratio, and root chord as well 
as stiffnesses, tuning masses, and tuning mass locations. 
Aerodynamics and dynamics constraints are enforced. 
Aerodynamic constraints consist of limits on power required 
(in hover, forward flight and maneuver), on stall, trim, and 
drag divergence Mach number. Dynamic constraints are on 
frequencies, minimum autorotational inertia, and maximum 
blade weight. 

The procedure is demonstrated for the design of a scale 
model of a rotor blade for a utility helicopter. Results are 
presented for two cases. Case 1 has an objective function 
involving performance (in hover, forward flight, and maneu- 
ver) and forward flight dynamics. Case 2 has an objective 
function involving performance in hover and forward flight 
dynamics. In the first case, the procedure is able to obtain a 
design which meets all the design requirements while reduc- 
ing the power required in hover and forward flight by 7 per- 
cent and in maneuver by 4 percent. The 4 per rev vertical 
hub shear in forward flight is reduced by 14 percent. This 
design compares favorably with predicted performance of the 
reference blade for the given flight conditions. In the second 
case, the optimized design improves the power required in 
hover and forward flight by 4 percent with no change in the 
power required in maneuver. However, there is a substantial 
reduction (68 percent compared to the initial blade) in for- 
ward flight hub shear at a cost of a 65 percent increase in 
maneuver hub shear. 

The designs from the integrated procedure are compared 
with designs from a sequential optimization approach in 
which the blade is first designed for performance and then 
optimized for dynamics. Using the Case 1 objective func- 
tion, the integrated approach is superior to tire sequential 
approach from the standpoint of achieving a better design for 
both performance and dynamics behavior. Using the Case 2 


objective function, the integrated approach design has 
somewhat worse performance measures. However, the inte- 
grated design has significantly lower 4 per rev forward flight 
vertical hub shear dynamic behavior than the sequential 
approach design. 
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Table 1. Design variable bounds, parameters, and flight 
conditions used in optimization examples 


F Design variables i 


Lower bound 

Upper bound 

Twist (deg) 

-20.0 

-5.0 

Taper initiation (r/R) 

0.265 

0.985 

Taper ratio 

0.05 

5.0 

Root chord (in) 

2.70 

8.10 

EI XX at root (lbm-ft 2 ) 

5000.00 

20000.00 

EI XX at y^- (lbm-ft 2 ) 

5000.00 

20000.00 

EI XX at tip (lbm-ft 2 ) 

5000.00 

20000.00 

EI ZZ at root (lbm-ft 2 ) 

5.00 j 

1000.00 

EI zz at yy (lbm-ft 2 ) 

5.00 

1000.00 

EI ZZ at tip (lbm-ft 2 ) 

5.00 

1000.00 

GJ at root (lbm-ft 2 ) 

5.00 

1000.00 

GJ at y tt (lbm-ft 2 ) 

5.00 

1000.00 

GJ at tip (lbm-ft 2 ) 

5.00 

1000.00 

mi (lbm) 

0.0 

0.50 

m2 (lbm) 

0.0 

0.50 

m3 (lbm) 

0.0 

0.50 

Y1 (r/R) 

0.24 

0.985 

72 (r/R) 

0.24 

0.985 

y3 (r/R) 

0.24 

0.985 

Parameters 

AT . 
A1 mm 

3456 lbm-in^ 

An 

0.1 

Ct min 

1.0 in 

Number of blades 

4 

Number of aerodynamic segments 
HOVT 

CAMRAD/JA 

19 

18 

Number of structural segments 

50 

Number of design variables 

19 

Power available 

20 hp 

Blade radius 

56.22 in 

Maximum blade mass 

3.5 lbm 

Af 

0.1 per rev 

Flight conditions 

Rotational velocity 

639.5 RPM 

(in Freon density of 0.006 
slug/ft 2 ) 

Hover tip Mach number 

0.628 


Hover 

Forward 

flight 

Maneuver 

C L 

0.00810 

0.00810 

0.00985 

C D 

- 

-0.000811 

-0.000596 

Advance ratio 

- 

0.35 

0.30 
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Table 2. Initial versus integrated optimization design - 
Case 1 objective function 



Initial 

Design 

Final 

Design 

Change 

(%)* 

Twist (deg) 

-9.0 

-16.0 


Taper initiation (r/R) 

— 

0.68 


Taper ratio 

1.0 

1.79 


Root chord (in) 

5.4 

4.38 


mj (lbm) 

0.0 

0.150 


m 2 (lbm) 

0.0 

0.126 


m 3 (lbm) 

0.0 

0.018 


yi (r/R 


0.396 


y 2 (r/R) 

“ 

0.395 


y 3 (r/R) 

- 

0.475 


Hover (hp) 

15.4 

14.4 

-6.5 

Forward flight (hp) 

13.5 

12.6 

-6.7 

Maneuver (hp) 

12.3 

11.8 

-4.1 

Weight "(lbm) 

3.05 

3.34 

+9.5 

BAKnWITyiiWBMi 

3411 

3561 

+4.4 

(per rev) 

2.87 

2.72 


f 4 (per rev) 

5.57 

5.28 


(per rev) 

8.57 

8.23 


f 6 (per rev) 

9.66 

9.21 


tj (per rev) 

17.28 

17.10 


^ (per rev) 

15.14 

15.79 


S 4 ff (lbf) 

2.19 

1.89 

-13.7 

S 4 m (lbf) 

0.95 

0.92 

-3.2 


Final - Initial 
Change = 


Table 4. Initial versus integrated optimization design - 
Case 2 objective function 



Initial 

Design 

Final 

Design 

Change 

(%)* 

Twist (deg) 

-9.0 

- 8.1 


Taper initiation (r/R) 

- 

0.80 


Taper ratio 

1.0 

1.36 


Root chord (in) 

5.40 

4.50 


mj (lbm) 

0.0 

0.113 


m 2 (lbm) 

0.0 

0.157 


m 3 (lbm) 

00 

0.039 


yi (r/R) 

- 

0.246 


y 2 (r/R) 

- 

0.294 



- 

0.414 


I Hover (hp) 

15.4 

14.8 

-3.9 


13.5 

12.9 

-4.4 


12.3 

12.2 

- 0.8 

1 Weight (lbm) 

3.05 

3.36 

+ 10.2 


3411 

3501 

+ 2 . 6 " 

f 3 (per rev) 

2.87 

2.79 


f 4 (per rev) 

5.57 

5.45 


fg (per rev) 

8.57 

8.28 


f 6 (per rev) 

9.66 

9.82 


ti (per rev) 

17.28 

17.10 


^ (per rev) 

15.14 

15.90 


s 4 ff (lbf) 

2.19 

0.71 

-67.6 

S 4 m (lbf) 

0.95 

1.57 

+65.3 


Final - Initial ,, 
Chan 8 e x 100 


Table 3. Comparison of integrated optimum design 
(Case 1) with reference blade 



Present Method 

Reference 

Design 

Twist (deg) 

-16.0 

-16.0 

Taper initiation (r/R) 

0.68 

0.8 

Taper ratio 

1.79 

3.0 

Root chord (in) 

4.38 

5.40 

Hover (hp) 

14.4 

14.8 

Forward flight (hp) 

12.6 

13.1 

Maneuver (hp) 

11.8 

11.8 

Forward flight hub 
shear (lbf) 

1.89 

1.52 
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Table 5. Comparison of integrated and sequential optimization design procedures 



INTEGRATED 

SEQUENTIAL APPROACH 
Step 1 Performance Step 2 Dynamics 

Objective function 

Case 1: eqn. 10 
Case 2: eqn. 11 

Case 1: eqn. 12 
Case 2: eqn. 14 

Case 1: eqn. 13 
Case 2: eqn. 15 

Design variables 

®tW’ ytr> c r^ c t> *T 
EI XX , EI ZZ , GJ @ root 

EIxx- EI ZZ , GJ @ y lr 
EI XX , EI ZZ , GJ @tip 
mj, m2, m3 

yi> y 2 > y3 

®tw» ytr> Cr^t* c r 

EI XX , EI ZZ , GJ @ root 

^xx’ ®zz» GJ @ y^j- 
EI XX , EI ZZ , GJ @tip 
mj, m 2 , m3 

yi> y* ys 

Constraints 

Aerodynamics 

Dynamics 

Aerodynamics 

Dynamics 


Table 6. Comparison of sequential versus integrated optimization designs using Case 1 objective function 


Final Design 


INTEGRATED 


Sequential 

Integrated 

Change (%)** 

Change (%)* 



Sequential Change = 


Sequential Final - Initial 


X 100 


* „ . . , Integrated Final - Initial ,, , nn 

Integrated Change = 6 — — X 100 
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Table 7. Comparison of sequential versus integrated optimization designs using Case 2 objective function 



Initial 

Design 

Final Design 

Sequential 
Change (%)** 

Integrated 
Change (%)* 

SEQUENTIAL APPROACH 
Step 1 Step 2 

Performance Dynamics 

INTEGRATED 

Twist (deg) 

-9.0 

-16.0 

-16.0 

-8.1 



Taper initiation 

(r/R) 

“ 

0.63 

0.63 

0.80 



Taper ratio 

1.0 

1.28 

1.28 

1.36 



Root chord (in) 

5.40 

4.04 

4.04 

4.50 



mi (lbm) 

0.0 

6.0 

0.072 

6.113 



m 2 (lbm) 

0.0 

0.0 

0.006 

0.157 




0.0 

0.0 

0.043 

0.039 



yi (r/R) 

- 


0.288 

0.246 



Y2 (r/R) 

- 

- 

0.353 

0.294 



■SlHi 

- 

- 

0.413 

0.414 



Hover (hp) 

15.4 

14.4 

14.4 

14.8 

-6.5 

-3.9 

Forward flight (hp) 


12.6 

12.6 

12.9 

-6.7 

-4.4 


12.3 

11.9 

11.9 

12.2 

i -3.3 

-0.8 

iPKni 

3.05 

3.05 

3.17 

3.36 

+3.9 

+10.2 

QHSMi 

3411 

3411 

3456 

3500 

+1.3 

+2.6 

f} (per rev) 

2.87 

2.82 

2.77 

2.79 



f 4 (per rev) 

5.57 

5.24 

5.14 

5.45 



fg (per rev) 

8.57 

8.58 

8.45 

8.28 



f 6 (per rev) 

9.66 

8.99 

8.90 

9.82 



tj (per rev) 

17.28 

17.30 

17.10 

17.10 



t 2 (per rev) 

15.14 

14.94 

14.90 

15.90 



S 4 ff (lb0 

2.19 

1.75 

1.75 

0.71 

-20.1 

• -67.6 

S 4 m (lbf) 

0.95 

1.39 

1.44 

1.57 

+51.6 

+65.3 


” Sequential Change = SequenU ^Hna l - In Maj x 1(x) 
* Integrated Change = Imegratej Final - Initial x 100 
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Fig. 1 Rotor blade design variables 


Fig. 3 Reference blade 
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Fig. 2 Optimization procedure 
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Fig. 4 Initial blade design 
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a) Chordwise bending stiffness distribution a) Chordwise bending stiffness distribution 




b) Fiapwise bending stiffness distribution 


b) Fiapwise bending stiffness distribution 
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c) Torsional stiffness distribution 

Fig. 5 Initial and final stiffness distributions 
for integrated optimization - Case 1 


c) Torsional stiffness distribution 

Fig. 6 Initial and final stiffness distributions 
for integrated optimization - Case 2 
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a) Chordwise bending stiffness distribution 
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b) Flapwise bending stiffness distribution 
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c) Torsional stiffness distribution 

Fig. 7) Stiffness distributions for integrated 
and sequential optimization - Case 1 


c) Torsional stiffness distribution 

Fig. 8 Stiffness distributions for integrated and 
sequential optimization - Case 2 
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